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Abstract 



Following a recent work by Yoshino, we study the aging dynamics of a di- 
rected polymer in random media, in 1 + 1 dimensions. Through temperature 
quench, and temperature cycling numerical experiments similar to the experi- 
ments on real spin glasses, we show that the observed behaviour is comparable 
to the one of a well known mean field spin glass model. The observation of 
various quantities (correlation function, "donation" overlap function) leads 
to an analysis of the phase space landscape. 
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I. INTRODUCTION 



The study of directed polymers in random media has triggered a lot of interest and of 
works (for a recent review see 0), since it is related to many fields, from the fluctuations of 
interfaces to quantum mechanical problems in a time-dependent random potential || , or 
the very topical problem of vortex lines in high temperature superconductors [|J . It has also 
connections with spin glasses, as has been shown by Derrida and Spohn ||, who studied a 
mean field version (on a Cayley tree) of the random polymer, showing the existence of a low 
temperature phase similar to the Random Energy Model of Derrida M. 

In finite dimensions (d + 1, with d transverse dimensions), the existence of a phase 
transition has been shown for d > 2 [0-|9| , whereas the system is always in a low temperature 
phase for d — 1, and has been called a "baby spin glass" in [ |IO|j . 

The dynamics of such model has not been much studied so far. A recent numerical 
work by Yoshino |lTJ has made clear the existence of aging for the directed polymer in 
random media in dimensions 1 + 1, with violation of time translation invariance and of 
the fluctuation-dissipation theorem. The observed behaviour for the correlation function 
is similar to the scaling properties for the simulations of a three dimensional spin glass 



model ||12|| . It is analyzed along the lines of a scenario similar to the droplets model for 
spin glasses [0,0 : the polymer moves in a network of "tubes" where its probability of 
presence (calculated at equilibrium with a transfer matrix method) is high. The network 
has a quite complicated spatial structure, and the tubes form loops of various sizes. The 
dynamics consists then of rapid fluctuations inside the tubes (acting as traps) combined with 
thermally activated jumps between different tubes. These thermal excitations are compared 
to the droplets excitations. Besides, the fact that the loops display a broad distribution of 
|T~5l] induces a broad distribution of relaxation times, and thus aging. 



sizes 



On the other hand, Cugliandolo, Kurchan and Le Doussal, following the study of the 
aging dynamics of a mean-field spin glass model |16j , and of a particle in an infinite dimen- 
sional random potential ||TT| , |T8f , have proposed an analytical treatment of the long times 
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off equilibrium dynamics of an elastic manifold embedded in an infinite dimensional space, 
in the presence of a quenched random potential |19[ (the statics of such a model has been 
studied by Mezard and Parisi with a replica variational Gaussian approximation, becoming 
exact in this infinite dimensional limit ||20| ). This corresponds here to an infinite d. The 
finite dimension of the manifold leads to the study of the relaxation of its Fourier modes 
k. The two-times correlation and response functions Ck{t,t r ) and rk(t,t') satisfy dynamical 
equations where two regimes can be separated, as for the case of the particle fl7lJl8ll ) : a sta- 
tionary regime, and a regime displaying aging, where the properties of equilibrium dynamics 
(namely, time translation invariance and fluctuation dissipation theorem) are violated. Var- 
ious equations can be written for the long times behaviour of the functions, for example, the 
measure of the violation of the fluctuation dissipation theorem X(C) does not depend on k. 

We will here focus on the same model as |TT| ; after describing this model, the used 
dynamics, we perform several numerical experiments for aging dynamics. Preparing the 
polymer in its ground state, we also show that it can display stationary dynamics, and 
compare the two kinds of dynamics in order to analyse the phase space landscape. 



A. Model 

The polymer is defined on a square lattice of linear size iV : it consists of N monomers 
lying on the sites {(i,Xi), i — 1, • • • , N}, and the Hamiltonian is : 

JV 

H[{xi}} = £ (\xi - a*_i| + V(i, Xi)) , (1) 

i=2 

where V is a Gaussian random potential, with zero mean and variance a, uncorrelated from 
site to site. The first term gives the elastic energy ; moreover, the steps |xj — are 
restricted to or 1, and one of the extremities is fixed, so that a Boltzmann measure can 
be defined : x\ = N/2. 

Transfer matrix method can be used to study the statics of this model |L0| . In particular, 
the Edwards-Anderson parameter can be calculated, following pi : we take two replicas 
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of the polymer, with the same realization of the disordered potential ; if Z(x,y,L) is the 

partition function for the pairs of polymers arriving at transverse coordinates x and y after 

L steps (or monomers), and if Z(x, y, q, L) is the partition function restricted to such pairs 

having an overlap q, it is possible to write recursion relations in L for Z and Y(x,y,L) = 

Y, q q Z(x, y, q, L), and therefore to evaluate 

1 J2 T „Y(x, y,L) , , 

lim - x ' y { . 2 

l^oc L Y, x>y Z(x,y,L) 

and to average it over disorder to obtain qEA- We will use this value to check some expected 
long time limits of dynamical quantities (see below). 

The directed polymer is evolving with Monte-Carlo dynamics in a heat bath in the 
following way : a monomer and a move are chosen at random, and the move is performed with 
probability min(l, exp(— (3AE)) (Metropolis algorithm), (3 being the inverse temperature, 
and AE the change of energy involved. One Monte-Carlo step consists in iV such tries. 

B. Numerical experiments 

The transfer matrix method allows us to find the ground state in a given realization of 
a potential ; we have studied the dynamics in two cases : the initial configuration is either 
taken at random, or as the ground state. The polymer, therefore coming from an infinite 
temperature or from a zero temperature thermalized state, is then free to evolve at the 
temperature of the heat bath. We then measure the evolution of the energy, the two-times 
correlation function defined by 



C{t w + 1, t w ) = (£2 8 Xi {t w +t),xi (t w )), (3) 



where (.) is a mean over thermal noise and the overline denotes a mean over realizations of 
the quenched disorder, and the overlap between two copies of the polymer evolving in the 
same realization of the potential, with the same thermal noise until t w , and then decoupled 
(a process called "donation" in [2~I|) ; if these copies are labeled by (1) et (2) , this 



overlap is 
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Q tw (t w + t,t w + t) = (E s x p {tw+ t),^(t w+ t))- ( 4 ) 

1=1 

Besides, we have performed temperature cycling experiments in the same way as for real 



spin glasses ||25|| . Most of the runs have been made with a polymer of length N = 500, and 
some with iV = 800 : no finite size effects were seen for the used simulation times. 



II. AGING DYNAMICS 
A. Quench at initial time 

The initial configuration of the polymer is chosen at random, with the constraint |x. 



= or 1 ; then it evolves at a fixed temperature. It has already been observed (TT 
that the correlation function C(t w + t,t w ) previously defined displays aging behaviour : it 
depends explicitly on t w and t (see figure ([[])) ; as the system ages, it becomes more rigid, 
in the sense that it evolves slower and gets away from himself always slower. 

For times t much lower than t w , the dynamics has the characteristics of equilibrium 
dynamics : C(t w + t,t w ) depends only on t (time translation invariance), and it has been 
checked numerically |TT| that the fluctuation-dissipation theorem is valid. Besides, we have 



also checked the validity of the relation Qt w (t w + t,t w + t) = C(t w + 2t,t w ) |2J]. In this 
regime, we have therefore a quasi-equilibrium dynamics. Some well-known spin glass models 
(like for example the p-spin spherical model [Jl6| , or the random manifold JT^JT9|] ) present a 



correlation function decaying from 1 to qEA, with a power law approach to the q E A plateau : 
q EA + At'". 

We have been able to fit the results of the simulations with such form, using the values 
of qEA obtained by transfer matrix method : this fits therefore use only 2 parameters, and 
not 3 (see figure ([|) for the values, and figure (|1|) for examples of the fits). The obtained 
values of v (typically in the range 0.1 — 0.2) are displayed in figure (|2|). 

For times t comparable or bigger than t w , the correlation function C(t w + 1, t w ) depends 
explicitly on t w : the dynamics is no more time translation invariant. For t ^> t w , it decays 



as a power law, 

C(t w + t,t w )~ f(t w ,T)(t/t w )-\ (5) 

We are therefore in presence of a weak-ergodicity breaking behaviour, lim^oo C(t w +t, t w ) = 
The obtained values of A are comparable to the values of v (see figure (Q)), and both 



exponents are increasing functions of temperature, like in real spin glasses . 



It is important to note that the v exponent is different from the x exponent studied by 



Yoshino [|11|] : this x is obtained by the scaling form 

C(t w + t,t w ) ~r x (6) 

for the t <^t w part, with a global form 

C(t w +t,t w )=t~ x $(t/t w ). (7) 

If we take the limit lim^^oo lim^oo we obtain the same behaviour, but the opposite order 
of limits, lim^oo lim^^oo C(t ?i , + t,t w ) yields 0, in contradiction with the expected static 
limit 

lim lim C(t w + t, t w ) = q E A, (8) 

t— too t m —>oo 

which is also obtained by the form q E A + At~ u . 

With numerical data it is however difficult to prefer one of these forms, and much longer 
simulations (much bigger values of t w ) would be necessary 

If we now look at the overlap between two replicas separated at t w , Qt w , it seems that 
this function has a finite limit at large times t (see figure (f|)), with a value compatible with 
the value of qEA, for big enough t w . This constatation puts this model in the class I of the 



classification of [24], which includes domain-growth models, and the p = 2 spherical p-spin 



model |23| : it indicates that the evolution in phase space takes place in "corridors" |2T 
of size qEA- On the contrary, the case of the manifold embedded in a infinite dimensional 
space yields lim^^oo lim^oo Qt w (t w + t,t w +t) = 0, and therefore belongs to type II, which 
probably indicates a much more complex phase-space landscape, and occurs also for example 
for the p-spin spherical model with p > 3 [0] . 
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B. Temperature cycling experiments 



A spin-glass quenched under its transition temperature, and then submitted to tempera- 



ture changes, shows a very puzzling behaviour (see [25] for a review, and references therein) 
that we briefly describe now. After a quench at time t = 0, the temperature cycle is as 
follows (see figure (^)) : the temperature is T from t — to t — ti, then T + AT from 
t — ti to t — t 2 , and again T after t — t 2 ; AT can be negative as well as positive. For a 
positive AT, one observes a reinitialization of the dynamics, with, e.g. for a thermoremanent 
magnetization, a relaxation after t 2 identical to the one obtained after a quench at t = 
and a waiting time t e ff = t w — t 2 . On the contrary, for a negative AT, the system keeps 
memory of its evolution, and its relaxation corresponds to an effective age between t w and 

t w ~ t 2 + ^1- 

We have performed these jumps numerically with T = 2, AT = 1, 2, —1, and ti = 500, 
t 2 = 1500. We have then monitored the evolution of the energy of the polymer, as well 
as Cj urnp (t w + t,t w ) with t w = 2000, comparing these quantities with the ones obtained for 
AT = (constant temperature). 

It is clear (figure (|)) that Cj ump (t w + t, t w ) corresponds to a certain C(t e ff + t,t e ff), t e ff 
being an effective age for the system, depending on t 2 — 1 1 and on AT ; t e jf is less than t w 
for a negative AT, which means that the dynamic has been slowed down by the time spent 
at a lower temperature ; here we estimate t e // = 1000 = t w — t 2 + t\ (for AT = —1, the 
data corresponding to Cj ump (2000 + 1, 2000) with t± = 500, t 2 = 1500 are superimposed on 
the curve C(1000 + t, 1000)). For smaller values of AT, or longer times t 2 — t±, t e ff can be 
bigger than t w — 1 2 + ti : the time spent at T + AT can contribute a little to the aging. For 
positive values of AT, t e ff is bigger than t w (in figure (H), t e // = 3000), showing that the 
time spent at T + AT has contributed to the approach to equilibrium more than the same 
time at T. 

The behaviour of the directed polymer is therefore symmetrical for positive or negative 
variations of temperature. No reinitialization of dynamics is found. This type of behaviour 
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is similar to the one observed for mean-field spin glasses p5|,p3[ , and thus very different 
from the one observed in real experiments on spin glasses. 

It should be remarked that real experiments deal with response function, whereas we 
are monitoring correlation functions. However, numerical simulations of a three dimensional 



Ising spin glass model have shown also for the correlation functions [27] a partial reinitializa- 
tion of the dynamics for positive AT, and asymmetric outcomes of numerical experiments 
with positive or negative temperature cycles. Nothing of this kind is found here. 

III. RELAXATION FROM THE GROUND STATE 

The transfer matrix method allows to find the ground state of the polymer, given a 
realization of the potential. We then let the polymer evolve at temperature T as before. 

The measure of the correlation function C(t w + t,t w ) for various t w and t shows in this 
case a simple t dependence : the system is time translation invariant (figure (|7J). Besides, 
C(t w +t, t w ) does not seem to go to zero for large t, or at least stays well above the correlation 
at same times, for a system with random initial condition, for the accessible times : no power 
law decay t~ x is found. It seems that ergodicity is really, and not simply weakly, broken. 

We have also studied the dynamics for an initial configuration of energy close to the 
ground state energy, but spatially well separated. In this case, we observe a similar be- 
haviour, with a time translation invariant correlation function. 

Such behaviour has been observed for example in the TAP states of the p-spin spherical 



model with p > 3 [28 



IV. ENERGY 

We have monitored the evolution of the energy density E a i of the polymer ; for a random 
initial condition the initial energy is high, so the initial behaviour is a fast decay, followed 
by a much slower evolution. 



S 



When the polymer is prepared in its ground state on the contrary, its energy Ef grows 
quickly because of the thermal bath, and then stays constant. 

We show in figure @ the difference between the energy densities for both situations, in 
logarithmic scales : this plot shows that the evolution is compatible with a power law decay 
of E a i towards Ef. All the aging dynamics take therefore place at higher energy densities 
than those of the low lying states, where the dynamics is stationary. 



V. DISCUSSION 

Whereas the previous analysis of the dynamics made by Yoshino was based on an analogy 
with the droplets model []13|,|14|1 , we focus here more on a phase space analysis. 

It must however be clear that, as for many numerical simulations, especially for glassy 
systems, the available time scales remain quite small, and that the results should therefore 
be considered as tendencies, indications of behaviour. They allow us, with these precautions, 
to present the following analysis. 

The observed behaviour is quite similar to the one found for the spherical p-spin model 
with p = 2 : the dynamic consists in a slow search of the ground state, with a slowly decaying 
energy. There exists many states with low energy, but the polymer is not able to find them, 
during its aging dynamics : it remains at higher energy density ; on the contrary, if it is put 
in one of these configurations, it stays trapped and has a stationary dynamics. 

The behaviour of the overlap of two copies of the system, having same configurations 
until a certain time, and then decoupled, shows besides that the evolution takes place in 
some kind of "gutters" in phase space. The results of temperature cycling experiments also 
show that the directed polymer in 1 + 1 dimensions is a much simpler system than real 
spin glasses. These two results are probably related, and also in agreement, concerning the 



relative simplicity of the phase space, with the fact [|TT|] that the response to a tilt field 
applied at the end of the polymer does not display any aging ; however it would be nice 
to measure the response to a spatially sinusoidal field, and the relaxation after cutting the 
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field, to check if such a relaxation has aging, and how it depends from the wavelength of 
the field. However, the sample to sample fluctuations of the response functions are very 
important, and such measurements are therefore very difficult. 

These behaviours are in fact intermediate between the p = 2 spherical p-spin model 
(analogous to domain growth |23|] ) and the p > 3 case which, despite being a mean field 
model, displays a much more complicated behaviour with long term memory and a very 
complex phase space |16 |. 

It would certainly be very interesting to study the directed polymer dynamics in higher 
dimensions : these new dimensions could provide a way to avoid energy barriers by going 



around them. The appearance of these entropy barriers p9|-p2| could yield new interesting 
effects and a richer dynamic in a more complex phase space. 

It is a pleasure to thank M. Mezard who initiated this work and R. Monasson for useful 
discussions and comments. 
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FIG. 1. (a) : from bottom to top, C(t w + t,t w ) versus t for t w = 50, 500, 5000, 100000, 
at T = 1, and the fit C M (t) = 0.65 + .25 i~ - 08 . (b) : C(t w + versus t for T = 3, and 

*«, = 2000, 50000, 100000, and C as {t) = 0.24 + 0.545 t~ a185 . In both cases, q EA (T = 1) = 0.65 
and qEA(T = 3) = 0.24 have been obtained by transfer matrix method. 
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FIG. 2. Exponents A (for T = 1, 2, 3, 4, 5) and v (for T = 1, 2, 3) versus temperature. 
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FIG. 3. Values of qEA obtained by transfer matrix method, versus temperature. 
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FIG. 5. Temperature cycles. 
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FIG. 6. Results of the temperature cycling experiments : The symbols show C(t w + t,t w ) 
at constant temperature T = 2 for t w = 1000 (squares), t w = 2000 (crosses), t w = 3000 (dia- 
monds). The temperature cycles are done with t\ = 500, t2 = 1500, t w = 2000, T = 2, and yield 
Cjump(t w + t, t w ) : the corresponding curves are superimposed on C(1000 + t, 1000) for AT = — 1, 
and on C(3000 + t, 3000) for AT = 1. 
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FIG. 7. C(t w + t,t w ) versus t, in logarithmic scales : comparison of the evolutions for various 
initial conditions. For each figure, the three lower curves correspond to random initial conditions 
((a), for T = 1, t w = 500, 5000, 100000, and (b), for T = 2, t w = 10000, 50000, 100000), while 
three curves corresponding to the system being prepared in its ground state at initial time are 
superimposed onto each other (with the same waiting times t w as for the lower curves), showing 
that C(t w + t, t w ) depends only on t in this case. 
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FIG. 8. (a) : evolution of the energies as a function of time, for T = 2, for random initial 
conditions (upper curve), or for a system prepared in its ground state (lower curve), (b) : E a i — Ef 
versus time, for T = 2, in logarithmic scale (symbols) ; the straight line corresponds to the power 
law t-- 15 . 
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